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We develop a systematic derivation for the Limber approximation to the angular cross-power 
spectrum of two random fields, as a series expansion in (£ + 1/2) -1 . This extended Limber approx- 
imation can be used to test the accuracy of the Limber approximation and to improve the rate of 
convergence at large £'s. We show that the error in ordinary Limber approximation is 0(£~ 2 ). We 
also provide a simple expression for the 2 nd order correction to the Limber formula, which improves 

QQ ■ the accuracy to 0(£~ 4 ). This correction can be especially useful for narrow redshift bins, or samples 

with small redshift overlap, for which the th order Limber formula has a large error. We also point 
out that using I instead of I + 1/2, as is often done in the literature, spoils the accuracy of the 

£NJ . approximation to 

<L> ■ 

m ; 

O 1 I. INTRODUCTION 

i— i Many observations in cosmology are observations of random fields (e.g. the cosmic microwave background (CMB) 

4^ anisotropy, fluctuations in the mass density or galaxy distribution, the weak lensing shear or convergence field, and 
21cm emission line fluctuations). A primary means to learn about the distribution and evolution of large-scale 

q , structure are through correlation functions of these fields, the simplest being the two point correlation function or 

%-[ ' its Fourier transform, the power spectrum. Many observations are given in terms of the angular correlation function 

C/3 ! wab (ji • n') or its spherical harmonic transform, the angular power spectrum Cab{&) 

C$ . 

wab (fi ■ n')={A(h)B(n')) = £ ^±lc AB {t)P t {fi ■ n') (1) 

T— I ■ t 7F 

■ where A and B are line-of-sight projections of the fields being correlated (e.g. the temperature anisotropy AT/T or 
the mass fluctuation dp/p), n, nf are unit vectors indicating the direction of observation and Pi are the Legendre 

, polynomials. 

If} • Calculations of angular power spectra give expressions in terms of several integrals which must be evaluated nu- 
| merically. The Limber approximation [Q and its generalization to Fourier space [g, |j| is a commonly used technique 

■ to simplify calculations. In implementing the Limber approximation one assumes small angular separations (or large 
multipolc moment £) and that some of the functions being integrated are more slowly varying than others. The Limber 

. approximation is powerful method to accurately estimate the magnitude and understand the analytic dependencies 
of the projected power spectra. Also, since the Limber approximation reduces the number of integrals numerical 
calculations are simpler. 

■ In this paper, we present a systematic derivation of the Limber approximation to the angular power spectrum as 
a series expansion in (£ + which is a rigorous generalization of a technique introduced in [ij. While the first 
term in the expansion is the usual Limber approximation, higher order terms can be considered as an extension. We 
apply this approximation to a few examples where keeping additional terms in the expansion might be desirable. 
The results presented here can be applied to the cross-correlation of two random fields whose Fourier space power 
spectra are isotropic. An analysis of the Limber approximation and a proposed alternative approximation for the real 
space correlation function is given in For another discussion of some issues related to the validity of the Limber 
approximation for lensing power spectra, the reader can refer to Appendix C of Q. 

In fjll], we present the derivation of the extended Limber approximation for the angu lar cross-power spectrum of 



two random fields. In §111, we make a comparison with the flat sky approximation. §IV applies the derived first 
and second terms in the Limber approximation to a few examples: the galaxy auto-power spectrum, the cross-power 
spectrum of two redshift bins with small overlap in redshift, and the cross-power spectrum of broad and narrow 
redshift distributions. Concluding remarks are given in 
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II. DERIVATION OF THE LIMBER APPROXIMATION 

We first develop the theoretical expectation value of the cross-correlation of two random fields, projected on the 
sky. Let us consider two random fields -A(x) and -B(x) with their Fourier transforms defined as 

A(k) = J d 3 x e- tk x A(x) and B(k) = J d 3 x e - lk x B(x). (2) 

These fields could be, for instance, the density fluctuation Sp(x)/p or the Newtonian potential <&(x). The cross- 
correlation power spectrum, PabQz) (which is assumed to be isotropic) is defined by 

(A(k 1 )B*(k 2 )) = (2n) 3 S 3 (k 1 - k 2 )P AB (kx). (3) 

The projections of A and B on the sky are defined using Fa and F B projection kernels 

A(h) = J dr F A (r)A(rh), and B(h) = J dr F B (r)B(rn). (4) 

Now, expanding A and B in terms of spherical harmonics, the angular cross-power spectrum, Cab{£) is defined as 

Cab(I) = (A em B* em ) 

/d 3 k 
j^Pab (fc) (4tt) 2 ^ (kn )j t (kr 2 )Y em (k) Y/ m (k) 

dridr2F A (ri)F B (r 2 ) / —^—jt(kri)ji(kr a )PAB(k) 

kdkP A B{k) I dr 1 f A (r 1 )J e+1 / 2 (kr 1 ) / dr 2 f A (r 2 )Je+i/2(kr 2 ), (5) 



where jg's are the spherical Bessel functions of rank £ and Yg^s are the spherical harmonics. In the last step, we have 
substituted the spherical Bessel functions in terms of the Bessel functions of the first kind, Jt+\/ 2 , and defined: 

F A (r) , . _ F B (r) 
Mr) = — -f^; f B {r) = — (6) 

At the next step, we will develop a series representation for the integral of an arbitrary function multiplied by the 
Bessel function. We will use the fact that the Bessel function (for v > 0) grows monotonically from zero at x = to 
x ~ v and starts oscillating rapidly afterwards, to write: 

r-OQ 

hm / e-<*-^f{x)J v {x)dx = B Q f{v)+B 1 f'{v)+B 2 f"{v)+Bzf"'{u) + ... (7) 
Jo 

Using the Taylor expansion of f(x) around x = v=l + 1/2, wc find: 

B n = — lim/ e-< x - v \x-v) n J v (x)dx (8) 



n 



i) 



K -^-Yimi— / e-< x - v ^J v {x)dx. (9) 
n! e^o de n 1 



The integral over the Bessel function is a standard Laplace transform, which has a closed form: 



<~ iW * - «- . , 4 + £ + * _ + 0(e , , (10) 



yielding: 



1 z/ 3 19i/ 

So = 1, Bt = 0, B 2 = --, B 3 = - - , B 4 = - , B 5 = — , 
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Therefore, we find 



Cab(£) 



dk k P AB {k) ( k- l f A {r) - —fl(r) 



vk 



G 



vk 



2 ' DV ' 6 

where kr = v = £ + f /2. Combining the parentheses and collecting terms of the same order in v one finds, 



[11) 



dr 



CabW = / -PAB(-)f A (r)f B (r) f - - 



f 



r 2 ffl(r) , /£(r)\ , r 3 //<[' (r) , /£'(r) 



/aW f B (r)J 6 V/A(r) / B (r) 



+ 0(^ 4 ) (12) 



After algebraic manipulations and some integrations by parts, the two parentheses in the 1/v 2 term can be combined 
to find: 



Cab(1) = J ^P A B{k)f A {r)f B {r)^l+ V -^- 



d\nf A d In f B 
d In r d In r 



s(k)-p(k) 



0(v~ 



where 



s(k) 



d\nP AB {k) 
din k 



p(k) 



k 2 [3P^ B (k) + kP^' B (k)} 
3P AB (k) 



(13) 



(14) 



Equations (|13j|14|) show the first systematic correction to the Limber approximation, which can be used to reduce 
the error in the approximation from £~ 2 to £~ A . Moreover, we can use relative magnitude of the sub-leading term 
in the expansion as a criterion for the convergence/reliability of the Limber approximation. We thus see that the 
convergence of the Limber expansion depends on both v = £ + 1/2 and the f A , f B . If the two kernels f A and f B are 
peaked at the same distance f, the 1/v 2 term is subdominant when v > f/max[ovi, a B ] where a A is the width of f A 
and a B is the width of f B . However, if f A and f B are peaked at different distances, say r A > r B + a B , where r A and 
r B are the locations of the maxima, truncating the expansion requires v > f(r A — r B )/a A a B E5[. 



III. FLAT SKY AND 1 + 1/2 

Let us now think about the 2D power spectrum in the flat sky limit (for a comparison with the angular power 
spectrum see also 0). To do this, we will use cartesian coordinates with xn the line-of-sight direction and xj_ the 
perpendicular direction, and integrate along the X|| direction 

A(xj.) = Jdx\\F A (x\\)A{x ± ,x\\), (15) 

and 



A(k ± ) = j d'xL±e~ lk ^ I dx {{ F A (x\\)A(x±,x\\), (16) 
so the 2D power spectrum will be given by 

(A(kx)B(kl)) = J dxp A {xf\) J dxF\F B (xft) j ^e ife " (x lW )p AB (k ± , fc||)(27r) 2 ^(k ± + k'J. (17) 



Expanding the power spectrum about fc|| = (and assuming P AB (k±, k\\) = P AB {^jk 2 _ + k 2 )) gives 
P (v \ I A w ( \w ( \t> (v Ji x 1 1 f d\nP AB dlnF A (x ll )d\nF B (x l] ) \ 



2X||fcj_ \ d\nk d\nx^ rflnx|| 

'(18) 

where the derivatives of P AB are evaluated at fen = 0. Notice that all of the P AB factors are independent of xn. How 
to compare this to the angular power spectrum? We expect 

£(£ + l)C e ^ k 2 ± P{k ± ) (19) 
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FIG. 1: Left Panel: The galaxy auto-correlation for a narrow redshift bin (a z = 0.01) at redshifts zq — 0.3 (upper black curves) 
and zo = 0.6 (lower cyan/gray curves). The solid lines are the exact Ce (Equation (|^)) the dotted lines are the th order 
Limber approximation, the dashed lines are the Limber approximation keeping the first 0(v~ 2 ) term in Equation (|l3|). Right 
panel: The difference between the Limber approximation at th and 2 nd order (in \jv) and the exact angular power spectrum 
in redshift bins at zq = 0.3 (black) and zo = 0.6 (cyan/gray). The radii of convergence for the Limber expansion are roughly 
at I ~ 15 and 30 respectively. 



for large £. Expanding -^P{v jr) in Equation ( |l3| ) about r = r where r is, for example, the peak distance of Fa(t)Fb(t) 
and comparing with Equation ( |l8| ) we can see that indeed 1(1 + l)Ce f» kj_P2D(k±) for £ + 1/2 = fk±. Notice that 
while for large ts, I + 1/2 ps at small Ps the factor of 1/2 actually makes a difference. Comparing the Laplacian in 
spherical coordinates with the Laplacian in Fourier space shows that indeed kf — ► \J £{£ + 1) = £ + 1/2 + 0(ljt) is 
the correct replacement. 



IV. EXAMPLES AND COMPARISON OF LIMBER AND EXACT RESULTS AT DIFFERENT ORDERS 

Here we consider a few examples of calculations of angular power spectra. For simplicity we will assume a spatially 
flat cosmology so that r = x( z ) is the comoving distance. In all plots we assume a ACDM universe with £l m = 0.27, 
D.A = 0.73, fib = 0.046 as the fractional densities of matter, cosmological constant and baryons, Hubble constant 
today Hq = 70 km/s/Mpc, scalar fluctuation amplitude as = 0.8, and scalar spectral index n s = 0.95. For the linear 
matter power spectrum we use the transfer function of H. 



A. Power spectrum of a narrow redshift bin 



From Equation ([L3|) we can see that keeping the first term in the Limber approximation is accurate only so long as 
the functions /^(r - ) and are slowly varying. Here we calculate the galaxy auto-power spectrum (for instance, 

U [To], |ll|). In calculating the galaxy auto-power spectrum (ignoring non-linear evolution), C gg (£) these kernels take 
the form 

f,{r) = *&™%&D(z) (20) 

where H(z) is the Hubble parameter, D(z) is the linear growth function with D(z = 0) = 1, c is the speed of light and 
W(z, Zq) is a normalized selection function centered at zq. If we assume that the linear galaxy bias 6=1, the power 
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FIG. 2: Left panel: angular cross-power spectra between two samples with Gaussian width a z = 0.05. Upper black curves are 
for the cross power spectrum between bins at zo = 0.3 and zo = 0.4, lower cyan/gray curves are for more widely separated bins 
with zo = 0.3 and zo = 0.5. Solid lines show the exact power spectrum, dotted the th order Limber formula and dashed the 
Limber approximation to 2 nd order in 1/v. Right panel: the difference between the curves shown on left. 



spectrum in Equation ( |l3| ) is just the mass power spectrum P(k). If the selection function is too rapidly varying, 
one will need to keep additional terms in the Limber approximation. To illustrate this we take W(z, z ) to be a 
Gaussian centered at zq with variance a], and compare the exact expression for C gg (£) , Equation (^|) with the Limber 
approximation to zeroth and second order in \ fl Equation (13). From the discussion in || we expect the expansion to 
diverge for v < v c ~ (r(za)/a r ). For a z = 0.01, this gives v c w 15 at zq = 0.3 and v c 30 at zo = 0.6 . Comparison 
of the exact Ci with the Limber expansion at different orders is shown in Fig. [l]. For a given width, er, the Limber 
approximation is clearly more accurate for smaller zq at a fixed £. Very roug hly, for I > 5r(z )/a r the th order 
Limber approximation is accurate to ~ 1%. 



B. Cross-correlation of populations with small redshift overlap 



Consider the cross-power spectrum between two source distributions with a small redshift overlap. Here, we will 
use two selection functions with the same width but different mean redshifts. Cross-correlating different redshift bins 
is a tool for calibrating photometric redshifts (see, for example [[l2], |l3|). Distributions centered at different redshifts 
are also present in galaxy-lensing cross-correlation (which would more accurately correspond to a very broad and 
a narrow redshift distribution; for a review sec Jl4|). We then use the expression given in Equation ( p0| ) for each 
sample, but allow the central redshifts zq to differ. Comparison of the Limber and exact calculation at different orders 
is shown in Figure |^. As we had argued in [H], we see that the Limber approximation is less accurate for more widely 
separated redshift bins. Consequently, in this case, including the 2 nd order correction in Equation ([l^) could lead to 
a significant improvement in the accuracy of the Limber approximation 



C. Cross-correlation of broad and narrow source distributions 



Here we consider the cross-power spectrum between two sources with different redshift distributions, for example a 
broad and a narrow source distribution. This is analogous to galaxy-lensing cross-correlation where the lensing weight 
function is broadly distributed and the galaxy selection function is narrow. The limit that one source distribution 
becomes extremely broad is also analogous to the galaxy-CMB cross-correlation. We use the same Gaussian selection 
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FIG. 3: Left panel: angular cross-power spectrum between two redshift bins centered at z = 0.4 with different widths cti = 0.10 
and <T2 = 0.01. Solid lines show the exact power spectrum, dotted the th order Limber calculation and dashed the Limber 
approximation to 2 nd order in 1/u. Right panel: the difference between the curves shown at left. 



functions from the previous sections, but allow the widths of the two distributions to differ. This calculation is shown 
in Figure ^|. Even though one redshift bin is narrow, since the other is broad the Limber approximation still works 
very well. 



V. CONCLUSIONS 



We have provided a series expansion for angular power spectra: The first term in this expansion gives the usual 
Limber approximation Jl], g], while the higher order terms are an extension to the approximation. The expression 
for the Limber approximation to second order in 1/v is given in Equation (|i~3|), higher order terms can be derived from 
Equations (^|), (|^) and ([H]). Figures [l] through || plot the accuracy of the Limber approximation at th and 2 nd order 
in l/£ for a few examples. The Limber approximation is less accurate for rapidly varying projection kernels /a and 
/b, or for /a and fs with small redshift overlap. The extended Limber approximation derived here can be applied to 
a variety of situations such as galaxy, weak lensing and CMB auto-correlations or to cross-correlations between the 
different projected distributions. 

It is also worth pointing out that, even in the th order Limber formula, replacing j/ = £+l/2by€ (as is often done 
in the literature), will increase the error from 0(£~ 2 ) to Therefore, simply using I + 1/2, as obtained in our 

systematic derivation, can significantly improve the accuracy of the approximation. 
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